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ABSTRACT 


The  of  segmental  motion  on  the  density-density  correlation 

function  of  n viscous  polymer  liquid  has  been  analyzed  using  a general- 
ized relaxation  equation  developed  by  Zwanzig  and  Mori.  It  is  shown  that  for 
polymer  liquids  of  high  viscosity  Brillouin  scattering  is  closely  asso- 
ciated with  the  structural  relaxation  associated  with  the  motion  of 
chain  segments.  A single  relaxation  time  theory  is  shown  to  yield 
good  agreement  with  the  experimental  results  on  polypropylene  glycol. 

The  torsional  motion  involving  a small  number  of  monomer  units  is  shown 
to  be  responsible  for  the  dispersion  and  attenuation  of  the  hypersonic 
wave.  The  fact  that  the  Brillouin  scattering  spectrum  of  a polymer  li- 
quid is  insensitive  to  the  change  of  molecular  weight  is  discussed. 

We  have  shown  that  temporal  modulation  of  the  spatial  second  moment  of 
the  intermolecular  or  inter-segmental  interaction  energy  is  responsible 
for  the  relaxation  process  involved  in  Brillouin  scattering. 


INTRODUCTION 

Measurements  of  viscoelastic  properties  of  undiluted  polymers  and 
polymer  solutions  can  provide  information  about  the  nature  and  rates 
of  the  configurational  rearrangements.  A knowledge  of  the  translation- 
al and  rotational  motion  of  polymer  molecules  and  their  segments  has 
been  important  in  elucidating  the  structure  of  macromolecules  in  the 
liquid  state. ^ 

Modern  developments  in  optical  technology  have  made  polarized  Ray- 
leigh-Brillouin  scattering  spectroscopy  one  of  the  most  useful  methods 
for  the  investigation  of  viscoelasticity  of  macromolecules  in  the  high 

9 

frequency  (10  Hz)  range. 

Recent  experimental  studies  have  shown  that  the  Rayleigh-Brillouin 
spectra  of  polymer  fluids  are  sensitive  to  the  molecular  configurational 
rearrangements  and  segmental  motion  of  flexible  polymer  chains.^  There 
is  at  present,  however,  no  procedure  or  guideline  with  which  quantita- 
tive information  about  the  polymer  dynamics  can  be  extracted  from  the 
experimental  data  using  this  technique.  The  field  of  inelastic  laser 
light  scattering  of  polymers  and  polymer  solutions  is  presently  in  a 
stage  of  development. 

In  the  low  frequency  range,  the  bead-spring  (Rouse-Zim)  model  pro- 
vides an  adequate  description  of  the  dynamic  properties  of  macromolecules 
in  solution.  At  high  frequencies  the  polymer  chain  dynamics  are  governed 
by  localized  motion  involving  only  a few  segments.  The  length  of  the 
polymer  chain  is  not  expected  to  greatly  affect  the  chain  dynamics.  In 
other  words,  one  does  not  expect  the  Brillouin  spectrum  of  polymers  to 
be  sensitive  to  the  change  of  molecular  weight.  This  was  indeed  observed 


experimental  ly.  ^ 
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However,  we  still  do  not  know  at  present  how  the  small  scale  mo- 
tion involving  very  few  segments  affects  the  Rayleigh-Brillouin  spec- 
trum of  a polymer  solution.  In  other  words,  we  do  not  know  how  to 
relate  the  hypersonic  data  of  a polymer  system  to  the  dynamics  of  seg- 
mental motion. 

To  analyze  the  type  of  information  that  one  may  learn  by  studying 
the  Brillouin-Rayleigh  spectrum  of  a polymer  fluid,  we  have  investi- 
gated theoretically  the  effect  of  polymer  segmental  motion  on  the  den- 


i 
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sity-density  time  correlation  function.  The  density-density  time  cor- 

« 

relation  function  is  proportional  to  the  Fourier  transform  of  the  ob- 
served Brillouin-Rayleigh  spectrum.  In  the  present  paper,  we  have 
analyzed  the  relaxation  process  involved  in  affecting  the  peak  fre- 
quency and  the  linewidth  of  the  Brillouin  spectrum  of  a polymer  li- 
quid. Finally,  we  have  used  the  linear  response  theory  introduced 
by  Zwanzig  to  analyze  the  motion  of  the  normal  mode  of  a bulk  polymer 
liquid.  We  have  interpreted  the  Brillouin  spectra  of  pure  polypro- 
pylene glycol  liquids  in  terms  of  this  new  theoretical  result. 

THE  DENSITY-DENSITY  CORRELATION  FUNCTION  OF  A I’OLYMER  FLUID 
Consider  a pure  polymer  fluid  which  is  made  up  of  a collection  of 
identical  polymer  segments.  Each  macromolecule  contains  n-segments 
and  the  illuminated  volume  contains  N molecules.  The  polarized  Bril- 
louin-Raylcigh  spectrum  is  determined  mainly  by  the  I'ourier  transform 
of  the  time  correlation  function  (-‘(t),  given  by'*, 


where  C{t)  is  the  qth  spatial  component  of  the  density-density  corre- 
lation function;  p and  p'  are  indices  denoting  the  macromolecules, 

/ and  j,  k denote  the  segments  in  a macromolecule.  Rj(t)  is  the  posi- 

tion (the  bead  position)  vector  of  segment  j and  polymer  p with  re- 
spect to  the  laboratory  coordinate  frame  at  time  t.  The  summation 
over  p,  p',  j and  k are  over  all  polymer  segments  inside  the  scatter- 
ing volume,  q is  the  wave  vector  of  the  excitation  mode  selected 


through  the  Bragg  condition  in  light  scattering.  The  amplitude  of  q 
is  equal  to 


where  is  the  wavelength  of  incident  light  in  vacuum,  n is  the  index 
of  refraction  and  0 is  the  scattering  angle. 

For  fluids  consisting  of  small  size  molecules,  one  can  analyze  C(t) 
in  terms  of  entropy  and  pressure  fluctuations.  These  fluctuations  can 
then  be  calculated  using  the  linearized  hydrodynamic  equations  with  fre- 
quency dependent  transport  coefficients  in  a manner  similar  to  that  of 
ordinary  low  viscosity  fluids  consisting  of  small  molecules.  As  the 
polymer  molecule  increases  in  size,  the  center  of  mass  motion  becomes 
very  slow  and  the  diffusive  entropy  fluctuations  then  dominate  the  time 
evolution  of  C(t).  Moreover,  in  ordinary  lov;  viscosity  fluids,  the  width 
of  the  Brillouin  doublet  is  predicted  to  increase  with  increasing  vis- 
cosity.® For  polymer  fluids  of  high  viscosity  this  prediction  would 
lead  to  such  wide  Brillouin  peaks  that  they  would  not  be  observable. 

In  Fig.  1 the  peak  frequency  and  linewidth  of  the  Brillouin  peak  of 
polypropylene  glycols  are  shown  as  a function  of  viscosity.  For  the 
experimental  data  at  two  temperatures  the  frequency  as  well  as  the  line- 


width  appear  to  be  independent  of  viscosity.  This  is  in  complete  dis- 
agreement with  theoretical  results.  In  fact  the  Brillouin  peaks  of 
highly  viscous  polymer  fluids  are  quite  sharp  near  The  discrep- 

ancy between  the  results  of  an  ordinary  fluid  and  a polymer  fluid  is 
believed  to  be  due  to  the  segmental  motion  that  is  present  in  the 
polymer  fluid.  The  segmental  motion  leads  to  the  breakdown  of  the 
Navier-Stokes  equation.  As  a result  the  usual  hydrodynamic  equations 
cannot  be  used  to  calculate  the  Brillouin  spectrum  of  a polymer  fluid. 

The  structural  relaxation  which  has  a dominant  effect  on  the  propagation 
and  absorption  of  the  hypersonic  wave  in  a liquid  consisting  of  macro- 
molecules, but  this  is  not  adequately  accounted  for  in  the  hydrodynamics 
equation. 

For  a polymer  fluid  consisting  of  flexible  polymer  segments,  the 
correlation  function  C(t)  is  expected  to  depend  strongly  on  tempera- 
ture, due  to  the  fact  that  the  structural  relaxation  time  of  a polymer 
fluid  depends  strongly  on  temperature.  At  high  temperature,  the  rate 
of  segmental  relaxation  is  comparable  with  the  frequency  of  the  hyper- 
sonic sound  wave;  in  this  case  the  segmental  motion  will  contribute 
significantly  to  the  linewidths  of  the  Brillouin  side  bands.  At  low 
temperature,  the  rate  of  segmental  fluctuations  is  slow  compared  with 
the  hypersonic  frequency,  but  fast  compared  to  the  thermal  diffusion 
rate.  In  this  case  the  segmental  motion  will  manifest  itself  as  a broad 
line  sui)crim|)osed  upon  the  central  Rayleigh  component  arising  from  the 
translational  diffusion.  Therefore,  the  time  evolution  of  C(t)  is  asso- 
ciated with  the  whole  region  of  tlio  Bri llouin-Rayleigh  spectrum. 

The  question  of  great  interest  is  what  are  the  effects  of  the  struc- 
tural relaxation  on  the  Bri llouin-Rayleigh  spectrum?  To  answer  this 
question  a microscopic  theory  on  structural  relaxation  is  in  order. 


In  order  to  develop  a microscopic  theory  for  the  Brillouin  scat- 
tering spectrum  of  a bulk  polymer,  one  needs  to  investigate  the  cor- 
relation function  C(t)  as  defined  in  l-q.  (1)  which  can  also  be  writ- 


Cft)  = <5p(q,t)6p  (q)s 

where  5p((i,t)  is  the  qth  spatial  Fourier  component  of  the  density  fluc- 
tuation an'd  is  given  by 

fiefq.t)  = Z r = z r = 


where  the  direction  of  q is  taken  to  be  parallel  to  the  Z-axis,  and 

each  polymer  segment  j of  polymer  p contributes  to  the  density  fluc- 
• -P 

tuation  by  a term  e ’ j.  To  calculate  C(t)  using  a microscopic  equa- 
tion, we  first  realize  that  the  density  fluctuation  does  not  occur  in- 
dependently. 'I'he  momentum  and  energy  density  fluctuations  couple  also 
to  the  density  fluctuation.  In  this  paper  we  suppress  the  coupling  of 
the  density  and  energy.  This  is  justified  because  our  primary  interest 
lies  in  the  study  of  the  Brillouin  doublet  in  which  the  effect  of  ther- 
mal conductivity  can  be  neglected  for  an  isothermal  system. 

Among  the  three  com])onents  of  momentum  density,  only  the  longitudinal 
component  in  the  direction  of  q is  coupled  to  Ap , as  the  coupling  to  the 
transverse  components  perpendicular  to  q is  prevented  by  symmetry.^  The 
expression  for  the  longitudinal  momentum  density  is  given  by 

g^(q,t)  = Z Z m 

" ' P j (5) 

where  m^^ . is  the  mass  of  j segment  of  polymer  p. 


-()- 


The  qth  mode  of  the  density  fluctuation  gives  rise  to  two  Brillouin 
peaks:  the  Stokes  and  anti-Stokes  peaks.  The  two  peaks  correspond  to 
acoustic  waves  propagating  in  two  opposite  directions.  We  show  in  Ap- 
pendix I that  the  two  acoustic  modes  with  opposite  frequency,  u , are 
represented  by  the  dynamic  variables  and  C which  are  given  by 


. X r-  * 


- /2  <|«p|2X  ■ <|g 


(6) 


Thus  the  density  fluctuation  6p  can  be  expressed  in  terms  of  and  C . 


6p  = — CC  + C ) 
/2  ^ 


(7) 


The  density-density  correlation  function  C(t)  is  then  given  by 


C(t)  = 2 <|6pl2>{<^^(t)?^  > + <C_(t)C.  > + > + 


<C  (t)e  >} 


(8) 


where  the  first  two  terms  on  the  right  hand  side  of  Eq.  (8)  corresponds 
to  the  two  acoustic  waves  propogating  in  opposite  directions.  The  last 
two  terms  vanish  at  t =0,  but  are  finite  at  t ^ 0 due  to  the  dynamic 
coupling  of  the  two  acoustic  modes  with  opposite  . These  can  contrib- 
ute to  the  intensity  to  the  central  Rayleigh  component.  These  can  be 
neglected  in  calculating  the  spectrum  of  the  Brillouin  doublet.  There- 
fore, calculation  of  the  Brillouin  spectrum  reduces  to  the  calculation 
of  the  correlation  functions  given  by 


(>(t)  = > 


(9) 


Since  both  and  C depend  on  the  same  q and  give  similar  physical  re- 
sults, in  1 
for  C^(t)  . 


suits,  in  the  following  we  shall  use  for  either  or  C , and  C^(t) 


-7- 


I 


Tn  principle,  if  one  has  a proper  expression  for  the  equation  of 
motion  of  C^^(t),  the  correlation  function  C^(l)  can  then  He  calculated 
according  to  the  method  of  irreversible  statistical  mechanics.  Formal- 
ly, the  equation  of  motion  for  ^(t)  is  given  by 


(10) 


where  L is  the  Liouville  operator  which  governs  the  time  evolution  of 
the  dynamic  variable  t^(t).  While  it  is  straiglit-forward  to  write 
down  a formal  expression  for  L in  terms  of  the  kinetic  and  potential 
energy  operators  in  a many  body  system,  in  practice  this  is  never  done 
in  a polymer  system  due  to  the  large  number  of  internal  degree  of  free- 
dom and  unknown  intra-  and  inter-molecular  potential  functions.  For  a 
dilute  polymer  solution,  Kirkwood's  generalized  diffusion  equation  or 
the  stochastic  Fokker-Planck  equation  has  been  employed  to  evaluate 
the  time  correlation  functions  associated  with  the  intrinsic  viscosity^ 
and  internal  viscosity®,  however,  an  appropriate  Liouvill  operator  for 
a bulk  polymer  has  not  yet  been  developed. 

Nevertheless,  even  in  the  absence  of  an  appropriate  Liouville 
operator,  useful  insight  on  the  behavior  of  the  linewidth  and  frequency 
shift  of  a Brillouin  spectrum  of  a polymer  fluid  can  be  obtained  from 
Eq.  (10)  using  the  linear  response  theory  developed  by  Zwanzig®. 

We  start  by  defining  a projection  operator  P such  that  when  it 
operates  on  a general  function  0,  it  yields 


I’C  - Z F,  'F  ( *>'’  -iir,  *> 

q q q q q 

where  (1  is  a function  of  t .md  may  be  time  dependent. 


(11) 


The  correlation  function  C^(t)  can  be  constructed  by  using  P ac- 


cording to 


PC„(t)  = i: 


’q' 


where  we  have  assumed  that  there  is  no  correlation  between  normal  modes 
of  different  wave  vectors. 

According  to  the  well-known  procedure^  , the  equation  of  motion 
for  can  be  shown  (using  h'qs.  (11}  and  (12)  to  follow  a general- 


ized Langevin  equation, 
9C(^(t) 


= i w ('  ft)  -/  dr  K(t)  C (t-T) 

q q / q 


where  the  frequenc)'  u ^ is  given  by 
oj  = <rb^  If;  *><r  r *>”' 

q ' q q q^q 

and  K(i)  is  the  memory  function  given  by 


where 


K(i)  = <r(T)F*(0)><^^^C^^*>"> 


1.^)  = ein-l’)LT 


l'(0}  = I-  = (l-P)  i b 5 


The  generalized  Langevin  equation  as  given  in  Hq.  (13)  is  very 
useful  because  of  its  mathematical  simiilicity  atul  case  of  application. 
This  etjuation  can  be  used  to  calculate  the  Brillouin  spectrum  when  the 
macroscopic  hydrodynamic  equations  fail.  Using  this  equation,  we  can 
translate  physical  ideas  about  the  nature  of  Brillouin  scattering  from 
a polymer  fluid  into  mathematical  terms  with  a minimum  set  of  assump- 


-0- 


tions  about  the  complex  nature  of  the  fluid  state  of  macromolecules. 

In  the  absence  of  an  explicit  Liouville  operator  to  calculate  w 

q 

and  K(t),  we  can  argue  that  for  a polymer  fluid,  there  are  strong 

interactions  between  translational  and  orientational  motions  of  each 
segment  and  between  segments  so  that  the  normal  modes  C should  equi- 
librate  rapidly.  Thus,  we  may  expect  a set  of  fully  quantized  vibra- 
tional states  with  an  energy  separation  approximately  equal  to  the  hy- 
personic frequency.  In  flexible  chain  molecules  at  ordinary  tempera- 
tures, the  vibrational  relaxation  rates  for  the  quantized  vibrational 
modes  are  fast.  Moreover,  there  will  also  be  a distribution  of  relax- 
ation times  for  chain  segments  to  change  from  one  state  to  another. 
Therefore,  for  the  normal  mode  given  in  fq.  (6),  the  frequency  u may 
take  a real  finite  value.  We  assume  that  is  equal  to  the  undamped 

hypersonic  frequency.  We  show  in  Appendix  I that  m is  related  to  the 

q 

isothermal  compressibility  of  the  polymer  liquid. 

The  memory  function  K(t)  decays  exponentially  with  a distribution 
of  relaxation  times. 


K(t) 


f(Tj  e 


-t/tj. 


dT 


(18) 


where  f(T  ) is  the  distribution  function  of  relaxation  times.  For  a 


single  relaxation  time  fCt^)  is  proportional  to  a delta  functi 


on 


6Ct^-Tr). 


Knowing  and  K(t)  , wo  can  now  solve  nq.(i:S)  by  the  Laplace 


transform  techni()ue, 
Define 


K(z) 


r 

J o 


C^jltl  dt 


K(t)  dt 


(19) 


(20) 


'3 
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I'quntion  (36)  can  now  be  used  to  analyze  the  experimentally  deter- 
mined Brillouin  spectra  of  bulk  polypropylene  glycol.  The  Brillouin 
spectra  of  this  polymer  li(iuid  at  various  molecular  weights  have  been 
accurately  measured  in  our  laboratory.  To  fit  the  experimental  spec- 
tra, one  note's  that  Hq.  3(i)  contains  tliree  adjustable  parameters; 

(the  unperturbed  acoustic  frequency),  A (the  amplitude  of  the  memory 
function  K(t)  at  t = o)  and  (the  relaxation  time  of  the  memory  func- 
tion). However,  the  nature  of  these  parameters  is  well  understood  and 
no  ambiguity  to  the  curve  fit  will  exist. 

It  is  easy  to  sliow  that  A corresponds  to  the  square  of  the  modula- 
tion amplitude  associated  with  the  torsional  motion  of  polymer  segments 
in  the  direction  perpendicular  to  the  propagation  direction  of  the  hy- 
personic wave. 

The  unperturbed  acoustic  freciuency  also  depends  weakly  on  tem- 
])eraturc.  According  to  13q.  (26),  the  function  has  a maximum  at 

w = u)  , which  is  related  to  w according  to 

P T 


P = '"q'i-v 


'■'q'’  * 


he  'U'eoiii'  eqii.ilit)  holds  because  ^ '"p'’ 1 ) ^ ' 1- 


1 


IaiikU  ion  (27)  g 
the  ilispersion  curve 
tion  of  temperature. 


i ves  the  di’uier'.ion  relation  for 

has  a reverse  S-shape  wlien  ui  is 

P 

(irovided  that  we  assume 


o-p.  For 
plotted  as  a 


< 1. 

func- 


(28) 

At  hi>>h  temperature  one  has  At_2  <<  1 and  u approaches  oj  . Thus  co 

R p " q q 

is  eipial  approximately  to  the  Brillouin  peak  frequency  at  high  tempera- 
ture when  the  amplitude  of  modulation  is  small  compared  with  the  relax- 
ation rate. 

IVe  have  used  the  least  s<|uarcs  fitting  procedure  to  fit  the  theo- 
retical peak  frequency  and  linewidth  as  computed  by  using  Hq.  (26)  to 
the  experimental  values  at  various  temperatures.  The  fit  was  carried 
out  by  first  generating  a spectrum  at  one  temperature  using  reasonable 
values  of  A,  and  1:^.  The  theoretical  spectrum  was  then  compared 
with  the  exjicrimental  one.  The  procedure  was  repeated  with  different 
values  of  A,  t and  F.  until  the  theoretical  spectrum  quite  resembles 
the  exjieri menta 1 one.  Next  the  spectra  were  generated  at  all  temper- 
atures, from  which  the  peak  frequencies  and  linewidths  were  read  and 
compared  with  the  experimental  values.  The  difference  between  theoret- 
ical and  experimental  values  is  then  minimized  with  a least  squares 
fitting  program  by  further  adjusting  the  values  of  A,  x , and  E . 

O 3 

In  I'ig.  2,  we  show  the  Brillouin  peak  frequencies  and  the  spectral 
linewidths  (the  half  width  at  half  height)  determined  by  the  curve  fit- 
tin*’  procedure,  also  including  the  experimental  results  between  280  and 
‘100  K for  comparison.  The  agreement  between  the  theoretical  and  ex- 
jierimental  results  is  quite  good  in  this  temperature  region,  consider- 
ing the  fact  that  only  three  fiarameters  (A,  E and  T ) were  used  to  fit 

3 0 
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two  separate  sets  of  data:  frequency  shifts  and  linewidths.  It  should 
also  be  pointed  out  the  agreement  between  the  experiment  and  theory, 
given  in  Figure  2 could  be  improved  if  we  allow  to  depend  on  tem- 
perature. As  shown  in  Appendix  I,  is  pro])ortinal  to  (px^p) ' , where  p and 

Xy  are  density  and  isothermal  compressibility  of  the  fluid,  respectively. 

Both  quantities  depend  on  temperature  slightly.  Increasing  temperature 
has  a tendency  of  decreasing  because  pXy  increases  slightly  with  in- 
creasing temperature.  As  a result,  the  peak  frequency  and  linewidth 
shown  in  Figure  2 will  be  depressed  downward  at  high  temperature  and 
pushed  upward  at  low  temperature,  thus  bringing  a better  agreement  with 
the  experimental  result  than  the  fit  using  a temperature  independent 
which  is  used  in  the  present  work.  However,  since  we  don't  have  complete  p and 
X.J,  data  for  PPH,  the  present  procedure  is  sufficient  to  demonstrate  the 
validity  of  the  theory. 


Below  280  K,  the  calculated  Brillouin  peak  frequencies  are  smaller 
than  the  experimental  values,  due  to  the  fact  that  the  relaxation  rate 
^ becomes  comparable  to  at  this  temperature.  In  this  situation, 
one  would  expect  strong  coupling  between  the  normal  mode  C and  the 
velocity  F,  and  we  can  no  longer  consider  F to  be  a fast  dynamic  vari- 
able to  modulate  both  must  be  considered  with  equal  weight.  More- 
over, since  is  about  equal  to  the  modulation  amplitude  »^,  in  this 
situation  the  precise  functional  form  of  the  memory  function  becomes 
important,  and  we  do  not  expect  that  F,q.  12-1')  is  valid  at  low  tempera- 


Wc  now  consider  in  more  detail  the  quantity  A,  which  is  a measure 
of  the  square  of  modulation  amplitude  associated  with  the  motion  of 
polymer  segments  in  the  direction  perpendicular  to  In  Apendix  II, 


In  Apendix  II, 
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*■ 

■A 


L 


we  show  that  A can  he  written  as 


A . ili  I E I I 

^ ""  p p-i  j SZ?9Z'.’ 

r I j j j 


<,2  - w 2} 


q 

(29) 


where  u is  the  total  potential  energy  of  the  system  and 


m is  the 


mass  of  monomer.  The  3kTq  /m  term  can  be  calculated  and  is  found  to  be 


only  about  one  cpiarter  of  the  value  of  io^2  for  ppf,  jo  the  temperature 


range  _ considered.  The  important  term  for  A thus  comes  from  the 
potential  energy  part. 


We  then  assume  that  the  range  of  intermolccular  interaction  po- 


1, 


tential  is  short  compared  with  the  wavelength  of  light  (or  q ),  and 


in  this  case  we  can  expand  the  phase  factor  e^*^^^i  " in  powers  of 


q.  The  first  non-vanishing  q dependent  term  is  found  to  be  equal  to^*' 

= 


32,1 


A = V if-  Z Z E E < fzl’  - Z^ 

2 '<■  p p'i  j szI’dzP  * J 

II  J j 


q2  1 

2 Nm 


3^u 


p p'i  j 3 (Z?  - zV)  ^ 


(z'.’  - Z^')"> 
1 1 


(.^0) 


According  to  F;q.  (2S),  A is  positive,  and  thus  A^  must  be  greater 


than  the  2 


term,  this  is  consistent  with  the  stal)ility  condition  con- 


dition for  a physical  system  for  which  we  always  have 


32u 


3(ZP  - Z*?  )^ 
1 J 


> 0. 


One  notes  that  A vanishes  when  Z*?  = Z^.  Thus  the  existence  of  A 


1’  1 .1  P 

depends  on  the  presence  of  intermolecul ar  (or  inter-segmental)  inter- 


actions. ibis  result  is  im))ortant  in  the  case  of  dilute  polymer  solu- 
tion, in  which  the  polymer  dynamics  deiiends  only  on  the  polymer-solvent 
and  intractiain  interactions. 

Since  the  i ntermolc'ciil ar  interaction  u is  short-range,  it  may  be 
expanded  in  tei'ms  of  power  series  of  ttic  int ermol eciil a r (or  inter-seg- 
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mental)  ilistancc  along  the  direction  c|, 


II  = u + 
o 


8u 


■ . . 3(Z?  - Z*?  ) 
P.P  .i,J  1 J 


9^u 


J 


^ ~ .9(Z?  - 
PP  IJ  1 J ^ 


fZ?  - Z?')'  * 
1 J 


(31) 


Comparing  with  Eq.  (31),  we  may  consider  in  Eq.  (.30)  as  the  spatial 
second  moment  of  the  intermolecul ar  potential  energy.  For  t 0 this 
spatial  second  moment  changes  because  of  molecular  motion  which  modu- 
lates the  intermolecular  distance  and  relative  orientations.  For  mac- 
romolecules, the  center  of  mass  motion  is  slow;  the  torsional  motion 
involving  only  relatively  small  numbers  of  segments  is  thus  most  effec- 
tive in  the  modulation  of  the  second  moment. 

Our  computer  curve-fitting  result  in  PPG  indicates  that  the  modu- 
lation amplitude  /a  is  insensitive  to  temperature  variation.  The  value 
of  /K  is  also  found  equal  to  throughout  the  280  - 400  K temperature 
range.  On  the  other  hand  the  relaxation  time  Tj^  is  found  strongly  tem- 
perature dependent.  Using  Eq.  (28),  we  have  found  = 0.94  x lO"^^ 
sec  and  the  activation  energy  E = 3.1  ± (0.02)  kcal/molc.  Within  the 
entire  temperature  range  we  have  further  found  that  /K  is  less  than 
unity.  This  validates  the  second  step  of  Eq.  (27),  thus  indicating 
that  the  Brillouin  scattering  spectra  of  PPG  between  280  and  400  K may 
be  described  in  terms  of  a stochastic  process  of  fast  modulation. 

As  pointed  out  above,  the  relaxation  time  cannot  be  associated  with 
the  center  of  mass  motion,  but  rather  with  the  motion  involving  a small 
number  of  segments.  The  rapid  motion  of  a chain  segment  from  one  posi- 
tion to  another  will  involve  a local  reorganization  of  the  structure  of 
the  polymer  liquid.  The  microscopic  picture  is  apparently  implied  in 
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the  phenomenological  theory  of  structural  relaxation  in  a viscous  liq- 
uid developed  by  Isakovich  and  Chaban.*^  In  reference  11,  the  diffu- 
sion of  holes  taking  place  in  the  viscoelastic  liquid  is  considered 
to  be  the  mechanism  for  the  structural  relaxation  process.  In  our 
present  microscopic  theory  is  the  measurement  of  relaxation  times 
associated  with  the  modulation  of  the  spatial  second  moment  of  the  po- 
tential energy.  Relaxation  occurs  because  of  fluctuations  of  segmental 
positions.  Thus,  the  present  microscopic  theory  serves  to  provide  the 
Isakovich  and  Chaban  phenomenological  theory  with  a more  satisfactory 
statistical  foundation. 

Finally,  it  should  be  pointed  out  that  despite  the  fact  that  we 
expect  in  general  a distribution  of  relaxation  times  for  a polymer  sys- 
tem, the  single  relaxation  time  model  seems  to  fit  satisfactorily  the 
experimental  results.  This  may  be  understood  from  Eq.  (30)  because 
only  the  modulation  of  the  short-range  intersegmental  interaction  po- 
tential can  make  a contribution  to  the  memory  function  K(t).  The  dis- 
persion and  relaxation  in  the  Brillouin  scattering  spectrum  of  a bulk 
polymer  liquid  involves  only  short-range  localized  segmental  motions. 
The  distribution  of  relaxation  times  due  to  the  motion  of  an  entire 
polymer  chain  will  not  enter  into  the  picture.  This  suggests  that  the 
Brillouin  scattering  spectrum  of  a bulk  polymer  liquid  should  be  inde- 
pendent of  molecular  weight.  This  is  clearly  in  agreement  with  the  ex- 
perimental result  on  polypropylene  glycol  >3  « ** 

.SUMMARY  AND  CONCLUSION 

We  have  analyzed  the  effect  of  segmental  motion  on  the  density- 
density  correlation  function  of  a polymer  fluid.  The  density-density 


correlation  function  is  proportional  to  the  Fourier  transform  of  the 
Bril louin-Rayleigh  spectrum  which  one  can  measure  experimentally. 

We  have  utilized  the  linear  response  theory  to  analyze  the  Bril- 
louin  spectra  in  terms  of  the  dynamics  of  polymer  chain  segments  in 
bulk  polymer  systems.  We  have  shown  that  for  polymer  fluids  of  high 
viscosity  Brillouin  scattering  is  due  to  normal  modes  associated  with 
the  motion  of  chain  segments.  When  the  modulation  amplitude  is  small 
compared  with  the  relaxation  rate,  we  have  shown  that  a single  relaxa- 
tion time  theory  will  provide  a satisfactory  description  of  the  hyper- 
sonic frequency  and  Brillouin  spectral  linewidth  data  as  a function  of 
temperature.  By  comparing  the  theoretical  expression  with  the  experi- 
mental result  we  have  also  obtained  the  activation  energy  associated 
with  the  relaxation  time  and  the  modulation  amplitude.  The  latter 
quantity  is  then  used  to  delineate  the  mechanism  involved  in  causing 
the  dispersion  and  attenuation  of  the  hypersonic  wave.  It  is  found 
that  the  localized  motion  which  causes  modulatioii  of  the  second  moment 
of  the  intermolecular  interaction  potential  is  responsible  for  the  dis 
persion  and  attenuation  of  the  hypersonic  wave.  We  have  thus  estab- 
lished the  fact  that  Brillouin  scattering  probes  only  short-ranged 
localized  motion.  We  have  also  explained  why  the  Brillouin  scatter- 
ing spectrum  of  a viscoelastic  polymer  liquid  does  not  depend  on  the 
molecular  weight. 

Another  significant  result  in  the  present  work  is  that  the  Bril- 
louin spectral  linewidth  is  related  to  the  relaxation  of  a correlation 
function  of  the  velocity  of  normal  modes  associated  with  the  chain 
segments  of  different  macromolecules.  It  is  thus  clear  from  the  mathe 
matical  point  of  view  that  the  maximum  position  in  the  Brillouin  line- 
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width  versus  temperature  curve  will  not  necessarily  correspond  to  maxi- 
mum point  in  the  NMR  or  versus  temperature  data  because  the 

three  techniques  are  sensitive  to  different  types  of  relaxation  pro- 
cesses. We  expect  in  general  that  Brillouin  scattering  will  provide 
supplementary  information  which  cannot  be  obtained  or  extrapolated  from 
other  relaxation  techniques. 
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APPENDIX  I 


The  frequency  oi  and  the  isothermal  compressibility. 

Q 

The  frequency  given  by  Eq.  (14)  can  be  either  positive  or 
negative  depending  on  whether  substitutes  for  or  ^ . In  either 
case,  the  derivation  is  similar.  We  give  only  the  derivation  asso- 
ciated with  Substituting  Eq.  (6)  into  Eq.  (14)  we  obtain 


u)  = i{<(L  r-) > 

^ 2 <|6p|2>^  <|g,|2>''^ 


+ <(L  ^^)  

<l5p|2>'^ 


(I.l) 


Since  the  Liouville  operator,  L,  is  Hermitian,  Eq.  (1,1)  can  be  re- 


written as 


, <(iL6p)g  > <g  (iL5p)  > 

0)  = ^(-i)  r r + (i)  i r) 

^ <|6p|2>-^<|gJ2>^  <|6p|2>'^<|gJ2>^ 


Using  the  definitions  for  6p  and  g^  (Eqs.  (4)  and  (5)),  we  obtain 


(i) 

(ii  = q n 

<i6p|2>^ 


According  to  the  principle  of  equipaitition  of  energy,  we  have 


i <|g  |2>''5  = (N— )"'* 


where  N is  the  total  number  of  segments  (light  scattercrs)  within  the 
scattering  volume,  V. 

In  the  limit  of  small  q,  it  is  true  that^ 

I 1 « V;  = f\7n2l-T\/  = /'M/.l.'r  V.  rr 


<|5p|2>H=  (Vp^kTXj)  = (NpkTXj) 


(1,5) 


(A. 1.2) 


I 


where  p is  the  number  of  segments  per  unit  volume.  Xj  is  the  iso- 
thermal compressibility.  Combining  Eqs.  (1,3),  (1,4),  and  (1,5),  we 
obtain  the  final  expression  for  u : 


10  = 
q 


(pmXy) 


(1,6) 


where  C.^,  = (pmx.j,)  is  the  expression  for  the  isothermal  sound  veloc- 
ity^. Equation  (1,6)  is  exactly  the  definition  of  isothermal  sound 
frequency.  Similarly,  for  ^ , the  result  is; 

= -qS  (1.7) 

where,  as  in  (1,6),  both  q and  C.j,  are  the  absolute  values. 
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APPKNDIX  II 


Reduction  of  the  Expression  of  Modulation  amplitude  A. 


Since  is  normalized,  using  the  result  of  Appendix  1,  we  reduce 


Eq.  (25)  to 


A = {<(iLC^)(iLC^)  > - 0)^2} 


(11. 1) 


Substituti,ng  Eq.  (6)  into  the  above  equation,  we  obtain  a new  expres- 
sion for  A: 


1 <(iLg  )(iLg  ) > 

A = ^ 5 . ^ 2} 


(11.2) 


Using  the  definition  of  g^  and  the  principle  of  equipaition  of  kinet- 
ic energy  for  each  segment,  we  now  expand  the  first  term  of  Eq,  (11,2) 
into 


<(iLg^)  (iLg^)  > 


<ig,r> 


1 


Nmkf  ^ 


Z1 


^^’’zi  ^zj  ^ ^ q^(>^T)2) 


1 3 


(11.3) 


where  the  double  summation  I Z of  Eq.  (6)  is  replaced  by  the  single 

summation  I (or  E)  of  the  segments  over  the  whole  scattering  volume; 
i j 

the  momentum  of  each  segment  in  Z direction  m^.  Zj  is  replaced  by  ' 
This  can  be  done  provided  that  each  segment  has  the  same  mass  m. 

We  understand  that  the  ensembly  average  <X>  of  any  physical  quantity  X 
IS  equivalent  to  multiplying  X by  the  distribution  function  — = — and 


then  integrating  the  product  over  the  whole  phase  space.  We  shall  as- 
sume that  the  Hamiltonian  H has  the  form 


(11.4) 


P.2 

I 


**  ^ 2tn.  “^^1  ’ ’'2  ’ ■ ■ ■ 

1 1 


The  ensembly  averaging  of  each  term  in  Eq.  (11,3)  can  be  carried  out 
as  follows. 


i/ 


^-B  H . 3u.2  .N  ,N 
e (^)2  d r d p 

1 


= - <-^  > 
8 3Z.2 

1 


(II. 5) 


where  to  obtain  tlie  last  step,  we  have  carried  out  integration  by  parts. 
It  is  easy  to  obtain. 


L <P.‘*>  = 3q2(kT)2 


(11,6) 


Likewise,  we  consider 


i/ 


P.p'  e'^^^^i  ■ = - I 

1 J ■’I 


''dZ.'"3Z. 
1 J 


d r d p 


Z 62 


/ ^3Z.9Z. 

J 1 j 


-B  H . iq(Z.  - Z.)  ,N  .N 
e ) e "•  1 d r d p 


. I ^ ^ iq (Z.  - Z . ) 

6 '‘3Z.3Z.  > 

1 J 


ZB2  j f^^*^3Z. 


-B  11  - . , 3 -B  II  , 2 -B  H 1 

e ) - iq(^  e ) - q-^  e ] 


„iq(Z.  - Z.)  jN  ,N 
e 1 d r d p 


+ M 


^ iq(Z.  - Z.) 

= j''>+M 

. ,,  1 32u  ..  iq(Z.  - Z.) 

where  M = - < (-g^-  - -^  ) e 1 j-'> 

i j 


(11. 7) 

(11. 8) 


In  order  to  obtain  liq.  (11,7)  inteprations  by  parts  have  been  used 
repetitively. 

('onsider  next  the  fourth  term  in  bq.  (11,3) 


iqkT  <(I\  - P.)  e ‘ 

= hi  f 

ez  / '‘3Z.  ■ az.^ 

J 1 J 

^ iSL.  f 

62Z  / ^3Z. 

- --pV  ^ " 


e-BH  ,iq(Z.  - Z.) 


-SH  _i_  „-BH  iq(Z.  - Z.)  .N  .N 
' ^ ^ ^ ^ ^ J d r d p 

I 


(11,9) 


Combining  bqs.  (11,3),  (11,4),  (11,5),  (11,6),  (11,7)  and  (11,9), 

and  substituting  the  result  into  Hq.  (11,2),  we  obtain  the  expression 
for  A as 


2 'T  : Nm  '3Z.3Z. 
^ J 1 J 


- m 2) 
m (] 


(11,10) 


which  is  equivalent  to  P.q . (29)  in  the  text. 


Figure  Captions 

W Figure  1.  The  Brillouin  peak  frequency  and  linewidth  as  a function 

'‘V  shear  viscosity.  The  solid  points  are  the  frequency 

anti  empty  ones  are  the  linewidth. 

< Figure  2.  Comparison  of  the  theoretical  and  experimental  Brillouin 

peak  frequency  and  linewidth  data  as  a function  of  temperature. 
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